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Abstract 

For a set W C L p (T d ) of multivariate periodic functions on the torus T d and a given function 
if G L p (T d ), we study the approximation in the L p (T d )-norm of functions / G W by arbitrary linear 
combinations of n translates of ip and measure the error in terms of the quantity 

M„(W, <p) p := sup{ inf{ ||/ - £ c W {- - y,)|| p : q G R,y, G T d } : / G W}, 

where N n := {1,2,..., n}. We prove upper bounds for the convergence rate of M n (Up (T rf ), K r> d) P 
where Up(T d ) is the unit ball in the Korobov space Kp(T d ) and n r ^ d is the Korobov function. To 
obtain the upper bounds, we construct approximation methods based on sparse Smolyak grids. The 
case M n (U2 (T d ), K r .d)2, t > 1/2 is especially important since K% (T d ) is a reproducing kernel Hilbert 
space, whose reproducing kernel is a translation kernel determined by the Korobov function K r ,d- 
We also provide lower bounds for the quantity M n {U^{T d )) 2 := inf{M„(£/ 2 r (T d ), 95)2 : <p G L 2 (T d )} 
which gives information about the best choice of <p. 

Keywords Korobov space; Translates of the Korobov function; Reproducing kernel Hilbert space; 
Sparse Smolyak grid. 

Mathematics Subject Classifications (2000) 41A46; 41A63; 42A99. 



1 Introduction 

The d-dimensional torus denoted by T d is the cross product of d copies of the interval [—it, tt] with 
the identification of the end points. When d = 1, we merely denote the d-torus by T. Complex- valued 
functions on T d are identified with functions on M. d which are 2ir periodic in each variable. We shall 
denote by L p (T d ), 1 < p < 00, the space of integrable functions on T d equipped with the norm 

\\f\\ p :=(2^)-^^J/(x)|^ x y /P . (1.1) 
* Corresponding author. Email: dinhzung@gmail.com. 
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For vectors x := (xi : I G No) and y := (yi : I G N<j) in T d we use (x, y) := X)zeN d xiyi for the inner 
product of x with y. Here, we use the notation Nd for the set {1,2,... , d} and later we will use for 
the set {0, 1, . . . , d — 1}. Given any complex- valued integrable function / on T d and any lattice vector 
j = (ji : I G Nd) 6 Z d , we let /(j) denote the j-th Fourier coefficient of / defined by 



/(j) := (2^ d f /(x) Xj (x)dx, 



where we define the exponential function Xj at x G T d to be Xj( x ) = e^j' x ). Frequently, we use the 
superscript notation M d to denote the cross product of a given set B. Also, for notational convenience 
we allow No and "Lq to stand for the empty set. 

The convolution of two complex- valued functions f\ and /2 on T d , denoted by /i * /2, is defined 
at x G T d by equation 

(/i*/ 2 )(x) :=(27r)- d [ /i(x)/ 2 (x-y)dy, 

whenever the integrand is in Li(T rf ). We are interested in approximations of functions from the 
Korobov space Kp(T d ) by arbitrary linear combinations of n arbitrary shifts of the Korobov function 
K r> d defined below. The case p = 2 and r > 1/2 is especially important, since K^iT 11 ) is a reproducing 
kernel Hilbert space. 

In order to formulate the setting for our problem, we establish some necessary definitions and 
notation. For a given r > and a lattice vector j := (ji : I G Nd) G 7L d we define the scalar Aj by the 
equation 

Aj := II \a 

where 

, ^Z\{0}, 

1 , otherwise. 



A/ := 



Definition 1.1 The Korobov function K r ^ is defined at x G T d by the equation 

«r,d(x) := J^Ar 1 Xj(x) 



and i/ie corresponding Korobov space is 

Kp(T d ) :={f:f = Kr,d*9, 9 ^ L p (T d )} 

with norm 

\\f\\K^(J d ) '■= \\9\\p- 

Remark 1.2 The univariate Korobov function k t ^\ shall always be denoted simply by K r and therefore 
K r> d has at x = (xi : I G Nd) the alternate tensor product representation 

K r,d(x) = Y[ K r(xi) (1.2) 
because, when j = (j'j : Z G N^) we Ziaue i/iai 

^(x) = ^ A j r i Xj ( X ) = e n = n -e^w- 
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B2m{t) = t ^' 2m (1 - K 2m (2TTt)). 



Remark 1.3 If r > 1, then X^jez d Aj 1 < oo and so K r ^ is continuous on T d . Therefore, we conclude 
that K r p (T d ) C C{T d ). For r > 1/2 we get that K r ,d G L 2 (T d ) and therefore we also conclude that 
K 2 (T d ) C C(T d ). These are the only cases which we shall consider in this paper. 

Remark 1.4 Since K r ,d(j) 7^ for any j G 7L d it readily follows that \\ ■ \\K r (T d ) * s a norm. Moreover, 
we point out that the univariate Korobov function is related to the one-periodic extension of Bernoulli 
polynomials. Specifically, if we denote the one-periodic extension of the Bernoulli polynomial as B n 
then for t G T, we have that 

2ml 

When p = 2 and r > 1/2 the kernel K defined at x and y in T d as X(x, y) := K2r,d( x — y) is 
the reproducing kernel for the Hilbert space K 2 (T d ). This means, for every function / G K r 2 {T d ) and 
x G T d , we have that 

/(x) =(f,K(-,x)) K r (Jd) , (1.3) 

where (•, -)i«;(T d ) denotes the inner product on the Hilbert space K 2 (T d ). For a definitive treatment 
of reproducing kernel, see, for example, pQ. 

Korobov spaces K p (T d ) are important for the study of smooth multivariate periodic functions. 
They are sometimes called periodic Sobolev spaces of dominating mixed smoothness and are useful 
for the study of multivariate approximation and integration, see, for example, the books [16] and |13j . 

The linear span of the set of functions {K 2 r,d(- — y) : y 6 is dense in the Hilbert space 

K 2 (T d ). In the language of Machine Learning, this means that the reproducing kernel for the Hilbert 
space is universal. The concept of universal reproducing kernel has significant statistical consequences 
in Machine Learning. In the paper [12J, a complete characterization of universal kernels is given in 
terms of its feature space representation. However, no information is provided about the degree of 
approximation. This unresolved question is the main motivation of this paper and we begin to address 
it in the context of the Korobov space K 2 (T d ). Specifically, we study approximations in the L 2 (T d ) 
norm of functions in K 2 (T d ) when r > 1/2 by linear combinations of n translates of the reproducing 
kernel, namely, K2 r ,d{- — yi), yi G T d , I G N n . We shall also study this problem in the space L p (T d ), 
1 < p < oo for r > 1, because the linear span of the set of functions {K r ,d(' — y) : y £ T rf }, is also 
dense in the Korobov space Kp(T d ). 

For our purpose in this paper, the following concept is essential. Let W C L p (T d ) and ip G L p (T d ) 
be a given function. We are interested in the approximation in L p (T rf )-norm of all functions / G W by 
arbitrary linear combinations of n translates of the function ip, that is, the functions (p(- — yz), yz G T d 
and measure the error in terms of the quantity 



M n (W,^) p := sup{inf{ ||/ - £ w (- - y,)[| p : q G R, yz G T d } : / G W}. 

The aim of the present paper is to investigate the convergence rate, when n — > oo, of M n (U p (T d ), K r ,d)p 
where U p (T d ) is the unit ball in ^'(T^). We shall also obtain a lower bound for the convergence rate 
as n — > oo of the quantity 

M n {U r 2 (T d )) 2 := inf{M n ([/£(T<V) 2 : <p G L 2 (T d )} 

which gives information about the best choice of ip. 
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The paper [10] is directly related to the questions we address in this paper, and we rely upon 
some results from [10] to obtain lower bound for the quantity of M n (Up (T d )) p . Related material can 
be found in the papers [S] and [UJ. Here, we shall provide upper bounds for M n (Up (T d ), n r ^)p for 
1 < p < oo,r > l,p ^ 2 and r > 1/2 forp = 2, as well as lower bounds for M^C/J (T d ))2- To obtain our 
upper bound, we construct approximation methods based on sparse Smolyak grids. Although these 
grids have a significantly smaller number of points than the corresponding tensor product grids, the 
error approximation remains the same. Smolyak grids and the related notion of hyperbolic cross are 
useful for high dimensional approximation problems, see, for example, [7j and [8]. For recent results 
on approximations and sampling on Smolyak grids see, for example, [3], [6], [H], and [15]. 

To describe the main results of our paper, we recall the following notation. Given two sequences 
{ai : / € N} and {&/ : I G N}, we write a/ <C h provided there is a positive constant c such that for all 
I G N, we have that ai < cb\. When we say that ai x b\ we mean that both ai <C hi and b[ <C a; hold. 
The main theorem of this paper is the following fact. 

Theorem 1.5 If 1 < p < oo, p ^ 2, r > 1 or p = 2, r > 1/2, then 

M n (U r p (T d ),^ d )p « n-^lognY^, (1.4) 

while for r > 1/2, we have that 

n- r {\ogn) r ^ « M n {U r 2 {T d )) 2 < ^(log™)^" 1 ). (1.5) 

This paper is organized in the following manner. In Section [21 we give the necessary background 
from Fourier analysis, construct methods for approximation of functions from the univariate Korobov 
space Kp(Y) by linear combinations of translates of the Korobov function K r and prove an upper 
bound for the approximation error. In Section [3j we extend the method of approximation developed 
in Section [2] to the multivariate case and provide an upper bound for the approximation error. Finally, 
in Section HI we provide the proof of the Theorem 11.51 



2 Univariate Approximation 

We begin this section by introducing the m-th Dirichlet function, denoted by D m , and defined at t G T 

as 

n u\ ST u\ sin((m+ l/2)t) 

D m (t) := ^ X l(t) = m (2.1) 

and corresponding m-th Fourier projection of / G L p (T), denoted by S m (f), and given as S m (f) := 
D m * f. The following lemma is a basic result. 

Lemma 2.1 Ifl<p< oo, then there exists a positive constant c such that for any m G N, / G KL~(T) 
and g G L p (T d ) we have that 

\\f-S m (f)\\ p < cm- r \\f\\ K r m (2.2) 

and 

\\S m (g)\\ p < c\\g\\ p . (2.3) 

Remark 2.2 The proof of inequality is easily verified while inequality 112. 3\) is given in Theorem 
1, page 137, of J^. 
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The main purpose of this section is to introduce a linear operator, denoted as Q m , which is constructed 
from the m-th Fourier projection and prescribed translate of the Korobov function K r , needed for the 
proof of Theorem Q3J Specifically, for / £ KUT) we define Q m (f), where / is represented as / = n r *g 
for g £ L P (T), to be 

Our main observation in this section is to establish that the operator Q m enjoys the same error 
bound which is valid for S m . We state this fact in the theorem below. 

Theorem 2.3 If 1 < p < oo and r > 1, then there is a positive constant c such that for all m £ N 
and f £ KI(T), we have that 

\\f-QmU)h < cm- r \\f\\ K r {T) 

and 

||Q m (/)|| P <c||/|| K , (T) . (2.5) 

The idea in the proof of Theorem 12.31 is to use Lemma 12.11 and study the function defined as F m := 
Qm(f) — S m (f). Clearly, the triangular inequality tells us that 

\\f-Qm(f)\\ P < \\f-S m (f)\\ p + \\F m \\ p . (2.6) 

Therefore, the proof of Theorem 12.31 hinges on obtaining an estimate for L p (T)-norm for the function 
F m . To this end, we recall some useful facts about trigonometric polynomials and Fourier series. 

We denote by T m the space of univariate trigonometric polynomials of degree at most m. That is, 
we have that T m '■= spanjxz : \l\ £ Z, m+ i}. The space T m is translation invariant, nested in the sense 
that C T s for m < s and is set theoretic multiplicative, that is, T m • T n Q T n +m- We also require a 
readily verified quadrature formula which says, for any / £ T s , that 



s \ s I 



Using these facts leads to a formula from [5] which we state in the next lemma. 

Lemma 2.4 If m,n,s £ N, such that m + n < s then for any fx £ T m and fi £ T n there holds the 
following identity 

'2ttA / 2nl\ 



i&s \ y \ / 

Lemma 12.41 is especially useful to us as it gives a convenient representation for the function F m . 
In fact, it readily follows, for f = n r * g, that 



where the function 9 m is defined — S m (K r ). The proof of formula (|2.8p may be based on 

the equation 

Sm( ^* 5) = 2^TT £ Sm{9) (^fl) {SmKr) (• " 2^Tl) ' (2 ' 9) 
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For the confirmation of (|2,9p we use the fact that S m is a projection onto Tm, so that 

S m {K r * g) = S m (K r ) * S m (g). 
Now, we use Lemma 12.41 with f\ = S m (g),f2 = S m (n r ) and s = 2m + 1 to confirm both (|2.8|) and 

(USD. 

The next step in our analysis makes use of equation (12.81) to get the desired upper bound for 
ll-^mllp- For this purpose, we need to appeal to two well-known facts attributed to Marcinkiewicz, see, 
for example, [18J. To describe these results, we introduce the following notation. For any subset A of 
Z and a vector a := (ai : I G A) and 1 < p < oo we define the / p (A)-norm of a by 

(Ei6A|Oil P ) 1/P } 1<P<00, 

sup{|a;| : / G A}, p = oo. 
Also, for each s G Z + , we introduce the set 

L s = {I:IgN,2 s <!< 2 s+1 - 1} (2.10) 
and the mapping W m : T m — > C 2m defined at / G T m as 



Lemma 2.5 If 1 < p < 00, then there exist positive constants c and d such that for any m G N and 
f G T m there hold the inequalities 

c\\f\\ p <(2m + ir 1 ^\\W m (f)\\ pMm+1 < c'\\f\\ p . 

Remark 2.6 Lemma \2.5\ appears in I18\j page 28, Volume II as Theorem 7.5. We also remark in the 
case that p = 2 the constants appearing in Lemma \2.5\ are both one. Indeed, we have for any f G T m 
the equation 

(2m + l)- 1 / 2 \\W m (f)\\ 2 ^ m+1 = ||/|| 2 . (2.11) 

Lemma 2.7 If I < p < 00 and there is a positive constant c such that for any vector a = (a,j : j G Z+) 
which satisfies for some positive constant A and any s G Z m , the condition 

and also ||a|| 00i 2 + < A, then for any functions f G L P (T), the function 

M a (f) ■■= E a ifc)Xi 

belongs to L p (T) and, moreover, we have that 

\\M a (f)\\ P < cA\\f\\ p . 
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Remark 2.8 Lemma 2.1 appears in JiffJ / page 232, Volume II as Theorem 4-14 an d ^ sometimes 
referred as the Marcinkiewicz multiplier theorem. 

We are now ready to prove Theorem 12.31 
Proof. For each j E Z we define 

bj := V S m (g) e "2^+T (2.12 

and observe from equation (|2.8j) that 

Fm = Yl b M~ r Xj- (2-13) 



Moreover, according to equation (|2.12p . we have for every j £ Z that 

b j+2 m+i = b j . (2.14) 
We also point out that the (2m + 1) x (2m + 1) matrix defined as 

1 / Inilj 

U = (e am+i : /, j G Z 2m+1 

is a unitary matrix. Therefore, we get for all I £ Z2 m +i that 

_) = 2 (2.15) 

and also using equation (|2.1ip applied to the function S m (g) and the fact that S m is a projection we 
get that 

J2 N 2 <Hsll!- (2-16) 
We split the sum on the right hand side of (|2.13p into two parts, namely, F m = + F~, where 

F m ■= Yl 3~ r bjXj^ F~ := Y -I ' h A r 

jeN\N m ieN\N ra 

Let us first consider the function F+. We decompose the set N\ N m as a disjoint union of finite sets 
each containing 2m + 1 integers. Specifically, for each j E N we define the set 

l mJ :={l:le N,j(2m + 1) - m < I < j(2m + 1) + m} 

and observe that N \ N m is a disjoint union of these sets. Therefore, using equation (|2.14p we can 
compute 

F m = Yj X/ ^ ^ Xl = ^2 Gm,jXj(2m+l)-m 

jen iei m ,j jeN 

where 

Gm, 3 := Y (l + 3^rn + l)-m))- r b lX i. 

2m+l 
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Hence, by the triangle inequality we conclude that 

ll-^mllp < y] \\G m ,j\\p- ( 2 - 17 ) 

Now, we shall use Lemmas 12.51 and [2 . 71 to estimate ||G m j|| p . For this purpose, we define for each j E N 
the components of a vector a = (on : I E Z_|_) as 



(Z + j(2m + l)-m))- r , IeZ 2m+1 , 
0, otherwise. 



and introduce the function G m := ^Cigz 2m+ i Since, we have that 

fi (l\ _ J ^3 ' ^ Z2 m +, 

Uml j ~ \ 0, otherwise. 

we may conclude that G m j = M a (G m ) and so Lemma [2~71 is applicable when a value of ^4 is specified. 
For a fixed value of j and m, we observe that the components of the vector a are decreasing and 
moreover, it is readily seen that oq < (jm)~ r . Therefore, we may choose A = (jm)~ r and apply 
Lemma 12.71 to conclude that 

||Gmj||p < p(jm)~ r \\G m \\p, (2.18) 

where p is a constant which is independent of j and m. Now, we apply Lemma [2.5l to the function G m 
(using the lower bound), observe by a straightforward computation that W m (G m ) = W m (S m (g)), once 
again, apply Lemma [231 (using the upper bound) to the function S m (g) and then appeal to inequality 
(]2.3p . Following these steps, we establish the existence of a positive constant r], which is independent of 
m, such that ||G m || p < ??||/||a'^(T)- We combine this inequality with inequalities (|2. 18|) and inequality 
(|2.17p to conclude, that there is positive constant 7, independent of m, such that H-f^Hp < 7||/||kj(T)- 
Likewise we may argue that H-^llp satisfies the same bound (perhaps with different constants). But, 
now, we recall that F m = + F~ and so inequality (12. 2p of Lemma 12.11 and inequality (I2.6P proves 
the first part of Theorem 12.31 

We now turn our attention to the proof of inequality (12.5p . Since n r is continuous on T, the proof 
of (|2.5p is transparent. Indeed, we successively use the Holder inequality, the upper bound in Lemma 
12.51 applied to the function S m (g) and the inequality (12. 3p to obtain the desired result. □ 

We end this section by providing an improvement of Theorem 12.31 when p = 2. 

Theorem 2.9 If r > 1/2, then there is a positive constant c such that for all m E N and f E K^T), 
we have that 

\\f-Q m (f)h < cm- r \\f\\ K r m . 

Proof. This proof parallels that given for Theorem 12.31 but, in fact, is simpler. From the definition of 
the function we conclude that 

Ml = E r 2r N 2 = E E fc ~ 2r N 2 - (2-i9) 

We now use equation (|2.14p to obtain that 

Ml! = E E (k + j(2m + l)-m))- 2r \b k \ 2 <m- 2r Y,r 2r E ^ ^ 
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and appealing to the inequality (|2.16p we finally get that 

Il^ + ll2 < IMll = ll/lll:;(T) 

which completes the proof. □ 



3 Multivariate Approximation 

In this section we extend the observations made earlier to the multivariate case. For our presentation 
we find it convenient to express the linear operator Q m defined in equation (|2.4p in an alternate form. 
Our preference here is to introduce a kernel H m on T 2 defined for x, t G T as 



1 ^ / 2iM / 2fl \ 



H m (x,t) - 2m+1 ^ ... v .,,„ + | 



and then observe when / = k t * g for g G L p (T) that 

Qm{f){x)= [ H m (x,t)g(t)dt. (3.2) 

The idea used in this section is to construct linear operators which act on multivariate functions 
fixing all variables but one and then applying Q m to the resulting univariate function. The precise 
description of what we have in mind starts with the following simple construction. For each I G N<f, 
x = (xj : j G Nd) G T d and t G T, we define a vector x;(i) G to have components given as 



xj, jGN d \{/} 
t, 3 = 1 



Thus, to form x;(t) we remove its Z-th component and insert in its place the variable t. So, of course, 
we have that Xi(xi) = x. Now, we define a collection ID; of complex- valued functions on T . A function 
/ is in provided that there is a w G L p (T d ) such that for every x G T d we have that 



/(x)= K r (xi - t)w(xi(t))dt. (3.3) 

Jt 

On the set H>i we define for each m G N a linear operator Qi m at / G Bj and x G T d as 

Ql, ro (/)(x)= / Hrnixi-^W^it^dt. (3.4) 
JT 

Alternatively, we may express Q/ m in the (perhaps, less convenient) form 

Qi,m(/)(x) = Q m («;(x I (-)))(xi). 

To get accustomed to our definition we shall demonstrate that the collection of linear operators 
Qm = {Qi,m '■ I G Nrf} are pairwise commutative on Kp(T d ). For simplicity of presentation, and 
without loss of generality, we concentrated only on the two operators Qi, mi and Q2,m 2 - We write a 
vector x G T d in concatenated form x = (u, v) where u = (xi, X2) G T 2 and v = (x2+j '■ j G N^_2) G 
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T d 2 . Also, we introduce a convenient subset Bi 2 of complex-valued function / for which there is a 
g G L p (T d ) such that x£T d 

/(x) = / «v, 2 (u- t)p(t,v)dt. 
Certainly i^(T d ) C Bi >2 C Di n D 2 and so we can compute Q 2 ,m 2 (f) at x G T d to be 

Q2,m 2 (/)( x ) = / H m2 (x 2 -t 2 ){ K r (xi-ti)g(ti,t 2 ,v)dti}dt2 
Jt Jt 

= / Kr(xi - tx){ / H m2 (x 2 - t 2 )g(ti,t 2 , v)dt 2 }dti. 
Jt Jt 

Therefore, we see that Q 2t m 2 (/) G ^1 and so we also have that 

(Ql,miQ2,m 2 )(/)( x ) = / H mi (xi -ti)H m2 (x 2 ,t 2 )g(t 1 ,t 2 ,v)dtidt 2 
Jt 2 

which implies that Qi >mi and Q 2t m 2 commute on Di ;2 . 

As was already pointed out in inequality (jl.2p . the multivariate Korobov function K T) d is a tensor 
product of d copies of its univariate version, and therefore, it readily follows that 

K;{Y d )<Z p|B,. (3.5) 

Consequently, we may represent the whole collection of linear operators Q m := {Qi t m '■ I G N^} 
simultaneously on i£T(T d ). For this purpose, we introduce a multivariate kernel Hi m on T 01 x T d 
defined for each x = (xj : j G N d ), t = (tj : j G N d ) G T d as 

Ify m (x,t) := H m (x h ti) } j K T (xj-tj). (3.6) 

ieN d \{z} 

Therefore, it readily follows for the function / = K r ^ * g where g G Lp(T d ) and x G T d that 

Qi,m(/)(x)= f ff/, m (x,t) 5 (t)dt. (3.7) 

Of course, when d = 1 Qi, mi reduces to Q mi . 

Our definition of the collection Q m of linear operators have been crafted so that they commute 
on the space Kp(T d ). A family of commuting operators allows us to use some simple results from 
Boolean Algebra. Specifically, for each lattice vector m = (m,- : j G N4) G N d we form the operator 

Qm • — Ql,mi- 

zeN d 

The operator Q m has a nice representation on i£T(T ). To describe it we introduce a kernel H m on 
T d x T d defined at x = (xj : j G N d ), t = (fy : j G N d ) G T d as 

tf m (x,t) := ^(x,t) (3.8) 
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and conclude for / G Kp(T d ) represented as / = k t ^ * g where g G L p (T d ) and x G T d we get that 

Qm(/)(X) = f fl- m (x,t)g(t)dt. 



Remark 3.1 iVoie £/iai t/ie multivariate kernel H[ m has two subscripts while the univariate kernel 
H m has only one subscript while the other multivariate kernel H m has a vector subscript. 

To assess the error in approximating the function / by the function Q m f we need a convenient 
representation for / — Q m f- Specifically, for each nonempty subset V C N<j and lattice vector m = 
(mi : I G Nd) G N d , we let |V| be the cardinality of V and define linear operators 

Q m ,V.= l[{I-Qi, mi ). (3.9) 

lev 

Consequently, it follows that 

T-Qm= Y, (- 1 )' Vhl Qm,V, (3.10) 



id 

where the sum is over all nonempty subsets of Nd- To make use of this formula, we need the following 
lemma. 

Lemma 3.2 If 1 < p < oo but p ^ 2 and r > 1 or r > 1/2 when p = 2, d is a positive integer and V 
is a nonempty subset o/N^, then there exists a positive constant c such that for any m = (rrij : j G N rf ) 
and f G Kp(T d ) we have that 

||Qm,v(/)[|p < rn ° w WfWm^y (3.n) 

vLLiev m i) 

Proof. First, we return to the univariate case and introduce a kernel W rm on T 2 defined at x, t G T 

W r , m (x,t) := K r (x -t)- H m (x,t). 
Consequently, we obtain for / = k t * g that 

f(x)-Q m (f)(x)= [ W r , m (x,t)g(t)dt. 

Therefore, by Theorem 12.31 and 12.91 there is a positive constant c such that the integral operator 
W r ,m '■ L p (T) — > L p (T) defined at g G L p and x G T to be 

Wr,m(jg)(x)= / W r , m (x,t)g(t)dt 
JT 

has an operator norm satisfying the inequality 

||W r , m || p := sup{||W r , m ((7)|| Lp : \\g\\ Lp < 1} < — r . (3.12) 

77V 

Our goal is to extend this circumstance to the multivariate operator Q mj v- We begin with the case 
that |V| = 1. For simplicity of presentation, and without loss of generality, we assume that V = {1}. 
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Also, we write vectors in M d in concatenated form. Thus, we have x = (xi,y),y G M. d 1 and also 
t = (ti, v), v G T d_1 . Now, whenever f = K r * g we may write it in the form 



f(x,y)= I K r (xi - ti)w(xx(ti))dti 

where 



T 



w(xi,y) := «r,d-i(y - v)p(xi, v)dv. 



By Theorem 2.3 and 2.9 we are assured that there is a positive constants ci such that for all y G T rf 1 
we have that 

y) - Ql,mi(/)(aJi,y)| p da:i < / M^i, y)| p cfei. (3.13) 

We must bound the integral appearing on the right hand side of the above inequality. However, for 
r > 1 we see that K r ,d-l G C(T d_1 ) while for r > 1/2 we can only see that G L2(T d_1 ). Hence, 

in either case, Holder inequality ensures there is a positive constant C2 such that for all x,t G T rf we 
have that 



Kzi,y)| p <ca / | 5 (xi,v)|Pdv. (3.14) 

We now integrate both sides of A3. 13f) over y G T^ -1 and use inequality (I3.14D to prove the result for 
|V|- 1. 

The remainder of the proof proceeds in a similar way. We illustrate the steps in the proof when 
V = N s and s is some positive integer in N^. This is essentially the general case by relabelling the 
elements of V when it has cardinality s. As above, we write vectors in concatenate form x = (y , y 2 ) 
and t = (v 1 , v 2 ) where y^v 1 G M s and v 2 ,y 2 G IR^ - *. We require a convenient representation for the 
function appearing in the left hand side of the inequality (I3.2p . This is accomplished for functions / 
having the integral representation at y\y 2 G T s , given as 



Z(y\y 2 )= / ^(y'-v'MvWv 1 , 



where g G L p (T d ). In that case, it is easily established for y 1 = (yf : I G N s ) and v 1 = (vj : I G N s ) 
that 

H(I-Qi,m i )(f)(y 1 ,y*)= I l[W r , m (yl,v})g(-v\y 2 )dv\ 

ieN s Jl ieN s 

This is the essential formula that yields the proof of the lemma. We merely use the operator inequality 
(|3.12|) and the method followed above for the case s = 1 to complete the proof. □ 

We draw two conclusions from this lemma. The first concerns an estimate for the L p (T rf )-norm of 
the function / — Q m (f). The proof of the next lemma follows directly from Lemma 13.31 and equation 

prop . 

Lemma 3.3 If 1 < p < oo but p ^ 2 and r > 1 or p = 2 and r > 1/2, then there exists a positive 
constant c such that for any f G Kp(T d ), m = (rrij : j G N d ) we have that 

\\f ~ Qm(f)\\ LPm < (min{m : C jGNd}) rll/lkW- 



12 



We now give another use of this Lemma 13.31 by introducing a scale of linear operators on KL (T d ) . 
First, we choose a k E Z + and define on Kp(T) the linear operator = I — Q 2 k. Also, we set T_i = /. 
Next, we choose a lattice vector k = (kj : j E N^) E Z d , set |k| = J2jeN d % an< ^ on K r p (T d ) we define 
the linear operator 

l€N d 

in the same fashion as we passed from Q m on Kp(T) to Q m on i*Tp(T rf ). 
The next lemma is an immediate consequence of Lemma 13.21 

Lemma 3.4 If 1 < p < oo, p ^ 2, r > 1 or p = 2, r > 1/2, £/ien i/iere is a positive constant c such 
that for any f E and k E there holds the inequality 

||T k (/)|| p < c2-'"l fe l||/||^ (Td) . 

In the next result, we provide a decomposition of the space Kp(T d ). To this end, for k E N we 
define on Kp(T) the linear operator Rk := Q 2 k — Q2^-i an d also set Rq = Q\. Note that for k E Z + 
we also have that Rk = T^-i — T^. Moreover, it readily follows that 

Q 2k = R i- 
i<=z k+1 

Let us now extend this setup to the multivariate case. For this purpose, if 1 = (lj : j E N^) and 
k = (kj : j E Nd) are lattice vectors in Z^_, we say that 1 -< k provided for each j E we have 
that lj < kj. Now, as above, for any lattice vector k = (kj : j E N^), we define on Kp(T d ) the linear 
operator i?k = rLeN d an d observe that 

Q 2k = ^i?i. (3.15) 

We now are ready to describe our decomposition of this space Kp(T d ). 

Theorem 3.5 If 1 < p < oo, p ^ 2, r > 1 or p = 2, r > 1/2, then there exists a positive constant c 
such that for every f E K p \T d ) and k E Z^ ; we have that 

\\RUf)\\ P < c2-'l k l||/||^ (Td) . (3.16) 

Moreover, every f E Kp(T d ) can be represented as 

f = fi k(/), ( 3 - 17 ) 

kez^ 

where the series converges in L p (T d )-norm. 

Proof. According to equation (|3.15p we have for any / E Kp(T d ) that 

Q 2 k(/) = I>i(/)- ( 3 - 18 ) 
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If each component of k goes to oo, then by Lemma [3731 we see that the left hand side of equation (|3.18j) 



goes to / in the L p (T d ) -norm. Therefore, to complete the proof, we need only confirm the inequality 



(|3.16p . To this end, for each nonempty subset V C and lattice vector k = (ki : I G N^) G N rf we 
define a new lattice vector ky £ N rf which has components given as 

Since = n«GN d an< ^ f° r ^ £ we have that Ri t ki = ^i,fc;-i ~~ ^J,Jfcj we obtain that 

^k = E (-l)' V| 7k v , 

VCN d 

and so by Lemma 13.41 there exists a positive constant c such that for every / G Kp(T d ) and k G 
we have that 

||iM/)llp< E ll*W(/)||p < c E 2- r|kv| ||/k ;(Td ) < c2-l k l||/||^ (ird) 

VCN d VCN d 

which completes the proof of this theorem. □ 

Our next observation concerns the linear operator P m defined on Kp(T d ) for each m G Z + as 

P m := £ i? fc . (3.19) 
|k|<m 

Theorem 3.6 If I < p < oo, p ^ 2, r > 1 or p = 2, r> 1/2, i/ien i/iere exists a positive constant c 
such that for every m G and f G -/^'(T^); 

||/-Pm(/)[|„ < c2- rm m d - 1 WfW^y 

Proof. From Theorem 13.51 we deduce that there exists a positive constant c (perhaps different from 
the constant appearing in the Theorem 13. 5j) such that for every / G Kp(T d ) and k G Z^, we have that 

\\f-Pm(f)\\p = II E R *(f)\\p < E H fl *(/)HP 
|k|>m |k|>m 

< c E 2- r|k| ||/||x ;( T d ) < c\\f\\ K r {Td) E 2- r|k| 

|k|>m |k|>m 

< c2-^m d - 1 ||/||^ (Td) . 



4 Convergence rate and optimality 

We choose a positive integer m G N, a lattice vector k G Zl with |k| < m and another lattice vector 

s = (sj : j G N d ) G %eN d Z 2*j +1 +i to defme tne vector yk,s = ( 2 fcj+i J +1 : J ^ N d ). The Smolyak grid 
on T d consists of all such vectors and is given as 

G d (m) := {y k>s : |k| < m,s G ®jm d ^ 2 k j +1 +l}- 
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A simple computation confirms, for m — > oo that 

\G d {m)\ = Y[(2 k]+1 + 1) ^2 d m d -\ 

|k|<mjeN d 

so, G d (m) is a sparse subset of a full grid of cardinality 2 dm . Moreover, by the definition of the linear 
operator P m given in equation f|3. 19|) we see that the range of P m is contained in the subspace 

span{K M (- - y) : y G G d {m)}. 

Now, we are ready to prove the next theorem, thereby establishing inequality (|1.4p . 

Theorem 4.1 If 1 < p < oo, p ^ 2, r > 1 or p = 2, r > 1 /2, then 

M n (U;,K r , d ) p « n-^lognr^ 1 ). 

Proof. If n G N and m is the largest positive integer such that |G d (m)| < n, then n x 2 m m d ~ 1 and by 
Theorem 13.61 we have that 

M n (U r p (T d ),K r4 ) < sup{||/ - P m (f)\\ Lpm : / G U r p (T d )} « 2^ m m d ~ 1 x n^logn)^- 1 ). 

□ 

Next, we prepare for the proof of the lower bound observation in Theorem 11.51 To this end, let 
P g (lR^) be the set of algebraic polynomials of total degree at most q onK 1 , and E m the subset of M. m 
of all vectors t = (tj : j G N m ) with components in absolute value one. That is, for every j £ N m we 
demand that \tj\ = 1. We choose a polynomial vector field p : R z — > W 71 such that each component 
of the vector field p is in P g (IR i ). Corresponding to this polynomial vector field, we introduce the 
polynomial manifold in R m defined as M m ^ i(? := p(lR'). That is, we have that 

M m ,^ := {(pj(u) : j G N m ) : Pj G ¥ q (R l ),j G N m , u G R 1 }. 

We denote the euclidean norm of a vector x in R m as ||x||2. For a proof of the following lemma see 

Lemma 4.2 (V. Maiorov) If m,l, q G N satisfy the inequality I \og( ie ™ q ) < ^, then there is a vector 
t G E m and a positive constant c such that 

inf{||t -x|| 2 : x G M mi ; i9 } > cm 1/2 . 

Remark 4.3 If we denote the euclidean distance oft G M m to the manifold M mj /„ by dist2{t,Wl m j s ), 
then the lemma of V. Maiorov above says that 

mp{dish(y,Mm,i, q ) : y G E ™} ^ cm ^- 
Theorem 4.4 If r > 1/2, then we have that 

n- r (logn) r ^ < M n {U r 2 ) 2 « n-^logn)^^ 1 ). (4.1) 
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Proof. The upper bound of (14. ip was proved in Theorem 14. 11 and so we only need to prove the lower 
bound by borrowing a technique used in the proof of \10\ Theorem 1.1]. For every positive number a 
we define a subset H(a) of lattice vectors given by 

M(a) := |k : k = (kj : j G N d ) G Z d , N < a}. 

Recall that, for a — > oo, we have that |H(a)| X a(loga) rf_1 , see, for example, [3J. To apply Lemma 
14.21 we choose for any n G N, q = ["(log n)~ d+2 \ + 1, m = 5(2d + 1) [n log re J and I = (2d + l)n. With 
these choices we observe that 

h q :=\M{q)\xm (4.2) 

and 

gxm(logm)^ +1 (4.3) 

as n — > oo. Also, we readily confirm that 

I ( 4emq\ 1 
hm — log — — = - 

n-+oo m \ I J 5 

and so the hypothesis of Lemma 14.21 is satisfied for n — >• oo. 

Now, there remains the task of specifying the polynomial manifold M mj z ;(? . To this end, we introduce 
the positive constant £ := q~ r m~ 1 / 2 and let Y be the set of trigonometric polynomials on T d , defined 
by 

Y: = {C E *kXk:t = (t k :k€lI(g))€E fc «}. 

keH(g) 

If / G Y is written in the form f = Q ^kee(g) *kXk 5 then / = K T) d*9 for some trigonometric polynomial 
g such that 

<c 2 £ |A k | 2 , 

k GH(g) 

where Ak was defined earlier before Definition ll.il Since 

C 2 £ lAkl^CVnH^m- 1 ^)!, 

keH(g) 

we see from equation (|4.2p that there is a positive constant c such that ||5||L 2 (T d ) — c f° r a ^ n £~N. So, 
we can either adjust functions in Y by dividing them by c or we can assume without loss of generality 
that c = 1. We choose the latter possibility so that Y C t/J (T"). 

We are now ready to obtain a lower bound for M^C/J (T d ))2- We choose any ip G L2(T d ) and let v 
be any function formed as a linear combination of n translates of the function <p. Thus, for some real 
constants Cj G K and vectors G T d ,j G N n we have that 

v = £ c i^(- - yj)- 

By the Bessel inequality we readily conclude for 

f = C *kXkGY 

heH(g) 



16 



that 



ii/-< (t ,)>c 2 E 

keM(q) 



c 



E 



(4.4) 



We now introduce a polynomial manifold so that we can use Lemma 14.21 to get a lower bound for 
the expressions on the left hand side of inequality (14. 4p . To this end, we define the vector c = (cj : j G 
N n ) G M n and for each j G N n , let Zj = (z^; : / G N^) be a vector in C d and then concatenate these 
vectors to form the vector z = (zj : j G N n ) G C nd . We employ the standard multivariate notation 



$ = n *5 



and require vectors w = (c, z) G x C nd and u = (c, Rez, Imz) G written in concatenate form. 
Now, we introduce for each k G 1HI(q) the polynomial defined at w as 



qfc(w) : = 



c 



E c i zk 



We only need to consider the real part of q& , namely, pk = Re q^ since we have that 



inf < 



/ k -^E 



c 



: G R,yj G T d ^ > inf <j ^ |t k - p k (u)| 2 :u£l' 

keH(g) 



Therefore, by Lemma l4.2l we conclude there is a vector t° = (t^ : k G G K hq and the correspond- 

ing function 

/° = c E 

keH( 9 ) 

for which there is a positive constant c such that for every v of the form v = YljeN n c iV 9 (" — Yj) we 
have that 

||/° - v\\ 2 V2(Jd) > <mh = q- r x n- r (logny( d -V 

which proves the result. □ 
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